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A  SIMPLE  MODEL  FOR  PREDICTING  THE 
FLUX-DISTRIBUTION  THROUGH  THE  FOCAL  PLANE  OF  A 
MULTIFACETED-CONCENTRATOR  SOLAR  FURNACE 

David  E.E.  Carlson,  Richard  B.  Diver,  and  Edward  A.  Fletcher 

Department  of  Mechanical  Engineering 
University  of  Minnesota 
111  Church  St.  S.E.,  Minneapolis,  MN  55455 

A  model  of  the  flux  distribution  at  the  focus  of  a  multi¬ 
faceted  concentrator  is  presented. 

We  have  described  a  small  doubly-reflecting  research  fur¬ 
nace. [1]  Its  optical  components,  a  heliostat  and  a  concentrating 
array  of  hexagonal -in-plan-form  mirrors,  focus  sunlight  at  the 
aperture  of  a  cavity  receiver.  A  perfect  paraboloid  of  revolution 
would  have  given  higher  concentration  ratios.  But  large  parabo¬ 
loids  are  difficult  to  make  and  manipulate.  Small  hexagons  are 
convenient  and  can  be  close-packed.  Spherical  mirrors  centered  on 
a  sphere  simplified  construction  and  were  economical.  In  this 
note  we  use  a  simple  model  to  describe  the  radial  distribution  of 
the  normal  flux  at  the  focus  of  such  an  array.  Related  problems 
have  been  considered  by  others.  (See  for  example  [2-5]). 

Figure  1  is  a  schematic  representation  of  a  cross-section  of 
the  concentrator  mirror  array  in  a  plane  which  includes  its  axis. 

In  our  idealization  the  sun  is  a  uniform  disc  which  subtends 
0.00931  radians.  A  solar  ray  bundle  reflected  from  a  concentrator 


mirror,  i,  to  the  focal  point  of  the  array,  F,  produces,  in  the 
plane  normal  to  it  a  circular  image  of  diameter 

(1)  «  2hisin2(<&i/4)  +  Bgf  , 

where  d>^  is  the  rim  angle  of  its  center,  h^  is  its  effective  dia¬ 
meter,  Bs  is  0.00931  radians,  and  f  is  the  focal  length  of  the 
mirror . 

In  the  plane  normal  to  the  array  axis  the  solar  image  from 
each  mirror  becomes  an  ellipse  of  diminished  intensity  whose 
minor  axis,  B^=D1#  is  also  given  by  equation  1.  Its  major  axis  is 
well  approximated  by 

(2)  Aa  =  D^cos^  . 

We  first  construct  a  flux  profile  through  the  plane  normal 
to  the  array  axis  assuming  perfectly  reflecting  surfaces.  The 
concentration  ratio  attributable  to  a  particular  concentrator 
mirror,  i,  is  given  by 

(3)  CA  *  Ani/Aei  -  h?cos(4)i/2)/4aibi 

where  Anj  ^  is  the  projected  area  of  mirror  i  in  the  plane  normal 
to  the  array  axis  and  Ae^  is  the  area  of  the  ellipse  in  the 
focal  plane  normal  to  the  array  axis,  ai=Ai/2,  and  biJ=Bi/2*Di/2. 

The  average  concentration  ratio  to  a  receiver  having  an 
aperture  of  radius  rj  is  the  sum  of  the  power  contributions  from 
each  of  the  mirrors  divided  by  the  unconcentrated  solar  power 
that  would  have  been  delivered  to  an  aperture  of  equal  area.  The 
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smallest  aperture  which  need  be  considered,  ^min^min'  is  the 
minor  axis  of  the  ellipse  corresponding  to  the  mirror  with  the 
smallest  rim  angle.  Within  this  radius,  the  flux  into  the  aper¬ 
ture  is  constant  and  at  its  highest  value. 


In  principle,  no  light  falls  outside  a  radius  equal  to  the 
major  axis  of  the  ellipse  from  a  mirror  having  the  largest  rim 
angle.  The  average  concentration  ratio  depends  on  the  aperture. 
It  is  given  by 


(4)  Cj  -  Pj/i»r?I 


where  I  is  the  solar  intensity.  With  n  facets,  the  power  through 
the  aperture  is 


(5) 


n 

l 

i=l 


1 3 


The  power  contributed  by  each  concentrator  mirror,  p^j,  is  given 
by 

(6)  Pij  =  FijIAicos(<)i/2)  , 


where  F^j  is  the  fraction  of  the  ellipse  from  mirror  i  captured 
by  the  aperture  of  radius  rj.  A^  is  the  cross-sectional  area  of  a 
concentrator  mirror.  For  convenience  we  divided  the  concentra¬ 
tor  mirrors  into  three  sets. 


Set  1:  rj  _<  b^.  If  the  aperture  is  small,  radius  rj  is 
included  within  the  ellipse  from  mirror  i.  Some  light  from  mirror 
i  falls  outside  the  aperture.  The  fraction  of  the  light  re¬ 
flected  from  mirror  i  which  enters  the  circle  of  radius  r  j  is 

,iV“n  by  ISM  I  Special 
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(’)  Fijx  «  T^/ail>i  . 

Set  2:  b^<rj<a^.  The  aperture  fits  within  the  major  axis  of 
the  ellipse,  but  it  extends  over  the  edges  of  the  minor  axis. 
Some  of  the  light  from  mirror  i  falls  outside  the  aperture,  and 
some  of  the  aperture  gets  no  light  from  mirror  i.  In  this  case 

(8)  Fij2  =  1  “  (2/7r)sin-1R  +  (2r?/iraibi)sin“1(biR/ri) 

where 

R  -  ((a?-r ? )/ (a^-b?))1/^ 

Set  3s  rj  >  a^.  All  the  light  from  mirror  i  enters  the 
aperture  and 


(9)  Fij3 


1  . 


The  total  power  into  the  aperture  having  the  radius  r j 
n  n 


(10)  Pj 


I  Pi  -  I  IFijAicos«t>i/2)  . 
i»l  i=l 


is 


Since  our  concentrator  has  a  six-fold  center  of  symmetry,  each 
rim  angle  has  either  a  six-  or  twelve  fold  degeneracy,  substan¬ 
tially  simplifying  the  summation  required  in  equation  10. 


An  example  solution  for  our  furnace  is  shown  in  Table  1.  The 
radius  of  the  aperture,  rj  is  1.8cm.  N  is  the  number  of  mirrors 
having  the  rim  angle,  its  degeneracy.  Ani  is  the  projected 

area  in  the  plane  normal  to  the  furnace  axis  of  the  mirrors 
having  the  rim  angle  0^.  This  area,  multiplied  by  the  solar  flux 
intensity,  I,  gives  the  solar  power  incident  on  each  mirror.  The 
major-  and  minor-axes  of  the  elliptical  images  are  a^  and  bx.  Fxj 


is  the  fraction  of  the  incident  power  captured  by  the  aperture. 
The  next  column  gives  the  power  supplied  to  the  aperture  by  a 
single  mirror  of  rim  angle  4)^  divided  by  the  power  incident  on 
the  normal  area  of  a  mirror;  it  is  the  effectiveness  with  which 
any  particular  mirror  delivers  power  to  a  3.6cm  aperture.  The 
last  column  gives  the  concentration  ratio  furnished  by  each  of 
the  mirror-sets,  and  its  sum  is  the  concentration  ratio  of  the 
furnace . 

Table  2  gives  the  variation  with  Tj,  of  the  average  concen¬ 
tration  ratio,  Cj,  in  the  region  bounded  by  rj,  of  an  optically 
ideal  furnace  having  the  configuration  of  our  furnace.  Table  2 
also  gives  the  average  concentration  ratio,  C(rj),  in  the  annular 
ring  between  rj  and  the  next  larger  rj.  The  variation  of  C(rj) 
with  r j  is  the  flux  distribution.  The  theoretical  peak  concen¬ 
tration  ratio  of  our  furnace  is  about  12500  suns.  The  theoreti¬ 
cal  peak  concentration  ratio  of  a  perfect  parabola  of  revolution 
having  the  same  rim  angle  would,  for  comparison,  be  about  23000. 
Our  measured  peak  concentration  ratio  was  7100.  The  difference 
between  theoretical  and  measured  values  is  attributed  to  a  mea¬ 
sured  reflectivity  of  0.68  which,  alone,  would  have  reduced 

the  concentration  ratio  to  8500,  and  dispersion  due  to  mirror 
imperfections  and  misalignment  which  was  responsible  for  the 
balance  of  the  loss. 
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TABLE  1.  CONTRIBUTIONS  OF  MIRROR  SETS 
TO  THE  POWER  DELIVERED  TO  AN  IDEAL  FURNACE 
OF  APERTURE  3.6CM 


rad 

Ani 

cm^ 

ai 

cm 

bi 

cm 

Fij 

Pi/lAi 

NCi 

.0780 

438.7 

1.356 

1.351 

1 

.999 

258.6 

.1356 

438.0 

1.382 

1.370 

1 

.998 

258.2 

.1566 

437.7 

1.396 

1.379 

1 

.997 

258.0 

.2080 

436.7 

1.437 

1.406 

1 

.995 

514.8 

.2365 

436.0 

1.466 

1.425 

1 

.993 

257.0 

.2740 

434.9 

1.510 

1.453 

1 

.991 

256.4 

.2852 

434.6 

1.524 

1.462 

1 

.990 

512.4 

.3176 

433.5 

1.570 

1.491 

1 

.987 

255.5 

.3471 

432.5 

1.617 

1.520 

1 

.985 

509.2 

.3656 

431.7 

1.648 

1.539 

1 

.983 

509.0 

.4009 

430.3 

1.715 

1.579 

1 

.980 

253.6 

.4175 

429.5 

1.749 

1.599 

1 

.978 

253.2 

.4259 

429.1 

1.767 

1.610 

1 

.977 

505.9 

.4494 

428.0 

1.821 

1.640 

.9965 

.971 

502.9 

.4873 

426.1 

1.915 

1.692 

.9605 

.932 

241.2 

.4946 

425.7 

1.934 

1.702 

.9512 

.922 

477.4 

.5091 

424.9 

1.974 

1.723 

.9310 

.901 

466.4 

.5374 

423.3 

2.057 

1.767 

.8860 

.854 

442.2 

.5712 

421.3 

2.165 

1.821 

.8211 

.788 

203.9 

.5779 

420.8 

2.188 

1.833 

.8076 

.774 

200.3 

.5780 

420.8 

2.188 

1.833 

.8074 

.734 

400.6 

.5978 

419.6 

2.258 

1.867 

.7681 

.672 

380.0 

.6302 

417.4 

2.382 

1.925 

.7063 

.625 

347.6 

.6555 

415.7 

2.487 

1.972 

.6602 

.603 

323.6 

.6681 

414.8 

2.543 

1.996 

.6380 

.592 

312.0 

.6742 

414.3 

2.570 

2.008 

.6275 

.560 

153.2 

.6929 

413.0 

2.658 

2.045 

.5958 

.560 

290.1 

.7301 

410.1 

2.847 

2.122 

.5361 

.501 

259.2 

.7423 

409.2 

2.914 

2.147 

.5175 

.482 

124.8 

.7482 

408.7 

2.948 

2.160 

.5086 

.473 

245.1 

.7669 

407.2 

3.057 

2.201 

.4813 

.446 

231.0 

.7791 

406.2 

3.132 

2.228 

.4641 

.429 

111.1 

.7971 

404.6 

3.247 

2.269 

.4395 

.405 

209.7 

C^  *  10525 


TABLE  2.  VARIATION  OF  CONCENTRATION  WITH  APERTURE 


rj 

C3 

C(r  j ) 

1.0 

12534 

12532 

1.1 

12532 

12532 

1.2 

12532 

12532 

1.3 

12532 

12151 

1.4 

12479 

10148 

1.5 

12179 

08275 

1.6 

11706 

06728 

1.7 

11138 

05465 

1.8 

10525 

04430 

1.9 

09900 

03254 

2.0 

09278 

02780 

2.1 

08674 

02179 

2.2 

08097 

01590 

2.3 

07544 

01195 

2.4 

07026 

00891 

2.5 

06545 

00633 

2.6 

06099 

00466 

2.7 

05689 

00355 

2.8 

05315 

00243 

2.9 

04971 

00140 

3.0 

04654 

00076 

3.1 

04364 

00034 

3.2 

04097 

00007 

FIGURE  CAPTIONS 


Figure  1.  Schematic  diagram  of  a  cross-section  of  the  concentra¬ 
tor-mirror  array  in  a  plane  which  includes  the  axis  of  the  fur¬ 
nace.  The  inset  shows  a  side  view  of  the  focal  plane,  Pf,  and  the 
associated  plan  view  of  the  elliptical  image  of  the  sun  produced 
on  the  focal  plane  by  one  of  the  mirrors.  The  dimensions  indi¬ 
cated  on  the  figure  are  defined  at  the  appropriate  places  in  the 


